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We measured the ground-state static electric-dipole polarizabilities of Cs, Rb, and K atoms 
using a three-nanograting Mach-Zehnder atom beam interferometer. Our measurements provide 
benchmark tests for atomic structure calculations and thus test the underlying theory used to 
interpret atomic parity non-conservation experiments. We measured acs = drreo x 59.39(9)A , 
apib = dTreo x 47.39(8) A^, and ok = 47reo x 42.93(7) A^. In atomic units, these measurements 
are acs = 401.2(7), ORb = 320.1(6), and q;k = 290.0(5). We report ratios of polarizabilities 
acs/oiRb = 1.2532(10), q:cs/qk = 1.3834(9), and ORb/oK = 1.1040(9) with smaller fractional un¬ 
certainty because the systematic errors for individual measurements are largely correlated. Since 
Cs atom beams have short de Broglie wavelengths, we developed measurement methods that do 
not require resolved atom diffraction. Specifically, we used phase choppers to measure atomic beam 
velocity distributions, and we used electric field gradients to give the atom interference pattern a 
phase shift that depends on atomic polarizability. 

PACS numbers: 32.10.Dk,03.75.Dg 


I. INTRODUCTION 

Measurements of static electric-dipole polarizabilities 
serve as benchmark tests for ab initio calculations of 
electric-dipole transition matrix elements. These calcn- 
lations require understanding quantum many-body sys¬ 
tems with relativistic corrections, and there are many 
different methods that attempt to calculate these matrix 
elements in a reasonable amount of computing time [I]. 
Testing these methods is important because these matrix 
elements are used to calculate several atomic properties, 
such as lifetimes, oscillator strengths, line strengths, van 
der Waals interaction potentials and associated cross sec¬ 
tions, Feshbach resonances, and photoassociation rates. 
Measuring alkali static polarizabilities as a means of test¬ 
ing atomic structure calculations has been of interest 
to the physics community since Stark’s pioneering mea¬ 
surements in 1934 [2]. Static polarizabilities have been 
measured using deflection [2-5], an E-H gradient balance 
[ 6 , 7], times-of-flight of an atomic fountain [ 8 ], and phase 
shifts in atomic and molecular interferometers [9-12]. 

We measured the static electric-dipole polarizabilities 
of K, Rb, and Cs atoms with 0.16% uncertainty using 
a Mach-Zehnder three-grating atom interferometer [13, 
14] with an electric held gradient interaction region. We 
used the same apparatus for all three elements, so we can 
also report polarizability ratios with 0.08% uncertainty 
because the sources of systematic uncertainty are largely 
correlated between our measurements of different atoms’ 
polarizabilities. 

We compare our measurements to ab initio calcula¬ 
tions of atomic polarizabilities and to polarizabilities de¬ 
duced from studies of atomic lifetimes, Feshbach reso¬ 
nances, and photoassociation specroscopy. We also use 
our measurements to report the Cs 6 pi /2 and 6 ^ 3/2 state 


lifetimes, Rb 5 pi /2 and 6 ^ 3/2 state lifetimes, and K 
and 4 p 3/2 state lifetimes and the associated principal 
electric dipole matrix elements, oscillator strengths, and 
line strengths. Then we use our measurements to report 
van der Waals Cg coefficients, and we combine our mea¬ 
surements with measurements of transition Stark shifts 
to report some excited state polarizabilities with better 
than 0.09% uncertainty. 

Testing Cs atomic structure calculations by measur¬ 
ing acs is valuable for atomic parity non-conservation 
(PNC) research, which places constraints on beyond-the- 
standard-model physics. The PNC amplitude Ernc due 
to Z°-mediated interactions between the Cs valence elec¬ 
tron and the neutrons in its nucleus can be written in 
terms of electric dipole transition matrix elements and 
the nuclear weak charge parameter Qw- Atomic struc¬ 
ture calculations are needed to deduce a value oi Qw 
from an Apnc measurement [15-18] to compare to the 
Qw predicted by the standard model [19, 20 ]. Our mea¬ 
surement of acs tests the methods used to calculate the 
relevant matrix elements and provides a benchmark for 
the ( 6 si /2 ||li|| Qpi/ 2 ) matrix element, one of the terms 
in the expression for Apnc- 

This is the Hrst time that atom interferometry 
measurements of polarizability have been reported 
with smaller fractional uncertainty than the pioneering 
sodium polarizability measurement by Ekstrom et al. in 
1995 [9]. This is also, to our knowledge, the first time 
atom interferometry has been used to measure Cs polar¬ 
izability. Because it is challenging to resolve Cs atomic 
diffraction—our nanogratings diffract our Cs atom beams 
with only 20 /xrad between diffraction orders—we de¬ 
signed an experiment with an electric field gradient in¬ 
stead of a septum electrode, such as was used in [9, 10]. 
We also developed phase choppers [22-25] to measure our 
atom beams’ velocity distributions instead of using atom 
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FIG. 1. (Color online) Diagram of the Mach-Zehnder atom interferometry apparatus with phase choppers (green) and electrodes 
that produce an electric field gradient (blue). Dimensions are shown in red and tabulated in Table I. The supersonic atom beam, 
shown in gray, is collimated by two slits si and s2 with widths wi and W 2 before entering the first grating. The nanogratings, 
gl, g2, and g3, are spaced longitudinally such that Li = L 2 , which causes an interference pattern to form at the position of 
g3. We consider the four separate interferometers in our data analysis, labeled with k = +2,+1,-1,—2, that form via 0th, 
±lst, and ±2nd order diffraction from gl. The atoms are detected by a platinum Langmuir-Taylor detector [21], indicated in 
orange. The pair of blue circles represents oppositely-charged cylindrical electrodes (extending perpendicular to the page) that 
form a virtual ground plane between them. The electric field from these electrodes polarizes the atoms and thereby shifts the 
interference pattern’s phase. The phase choppers are shown in green; each phase chopper is a charged wire next to a grounded 
plane. The geometry terms relevant to the pillars and phase choppers are displayed in Fig. 2 and discussed in Section II A. Due 
to the rotation of the Earth, the lab has a rotation rate about the vertical axis of fliab,?; = 38.88 /rrad/s that is also relevant to 
our analysis. 


diffraction to study velocity distributions, as was done 
in [9, 11]. These two innovations enable us to measure 
polarizabilities of heavy atoms such as Cs without re¬ 
solving diffraction patterns. Without the need to resolve 
diffraction, we can use larger collimating slits and a wider 
detector to obtain data more quickly. These innovations 
also reduce some systematic errors that are related to 
beam alignment imperfections. 

We improved the accuracy of our measurements com¬ 
pared to our previous work [11] by redesigning the elec¬ 
trodes that apply phase shifts to our interferometer. The 
new configuration of electrodes, two parallel, oppositely- 
charged cylindrical pillars, allows us to determine the dis¬ 
tance between the atom beam and the virtual ground 
plane between the pillars with reduced statistical uncer¬ 
tainty. We reduced systematic error by making more 
accurate measurements of the width of the gap between 
the pillars, the pillars’ radii, the voltages on the pillars, 
and the distance between the pillars and the first diffrac¬ 
tion grating. Our measurements also required a sophis¬ 
ticated model of the apparatus, which included interfer¬ 
ence formed by the 0th, ±lst, and ±2nd diffraction or¬ 
ders, the finite thickness and divergence of the beam, and 
the finite width of the detector [25]. Because beams of Cs, 
Rb, and K had different velocity distributions and diffrac¬ 


tion angles, we developed a more detailed error analysis 
in order to understand how those attributes affected the 
systematic uncertainties in polarizability measurements 
of different atoms. To support our error analysis, we also 
developed a method to monitor and adjust the distances 
between nanogratings in our interferometer. 


II. APPARATUS DESCRIPTION AND ERROR 
ANALYSIS 

A schematic diagram of the three-grating Mach- 
Zehnder atom beam interferometer we use to make our 
measurements is shown in Fig. 1. A mixture of He and Ar 
gas carries Cs, Rb, or K vapor through a 50 /rm nozzle to 
generate a supersonic atom beam [26, 27]. We adjust the 
carrier gas composition to change the beam’s average ve¬ 
locity: a higher percentage of Ar results in a slower beam. 
The atom beam passes through two collimating slits and 
diffracts through three silicon nitride nanogratings, each 
with period dg = 99.90(0.5) nm [28, 29]. The first two 
gratings manipulate the atoms’ de Broglie waves to form 
a 99.90 nm period interference pattern at the position of 
the third grating. The method of observing interference 
fringes is described in detail in [30]: we scan the second 
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TABLE I. List of apparatus dimensions described in Fig. 1 
and Fig. 2. Dimensions with no quoted uncertainty have un¬ 
certainty much less than what would be significant to our 
analysis, fld and 0^2 are the closest distances between the 
wires and the ground planes for phase choppers 1 and 2, and 
Upiiiars is half the width of the gap between the pillars (the 
closest distance between the virtual ground plane and either 
pillar). Li — L 2 = 0 ± 30 /rm, and the uncertainty in Li + L 2 
is insignificant. 


Zsl,s2 

860 mm 

'^gl,pillars 

833.5 ± 0.25 mm 

'S'gl,cl 

269.7 mm 

^s2,gl 

100 mm 

'^g2.c2 

598 mm 

Zcl,c2 

1269.3 ± 0.25 mm 

Ll 

940 mm 

L 2 

940 mm 

Wi 

30 ± 6 ^m 

W2 

40 ± 6 ^m 

Wdet 

100 ± 3 /rm 

^pillars 

1999.85 ± 0.5 Aim 

-^pillars 

6350 ± 0.5 Aim 

del 

986 ± 25 Atm 

Rci 

785.5 Atm 

ac2 

893 ± 25 Atm 

Rc2 

785.5 Atm 


Not to scale 



FIG. 2. (Color online) Diagram showing the dimensions, 
shown in red, that describe the static electric fields created by 
the pillars and by the phase choppers. The circle represents 
the cross-section of a metal pillar or charged wire with radius 
R held at voltage V. The GND line represents the ground 
plane, which may by physical (in the case of the choppers) or 
virtual (in the case of the pillars), a is the closest distance be¬ 
tween the pillar edge and the ground plane. The parameter b 
is the distance between the ground plane and the effective line 
charge within the pillar. The atom beam center is a distance 
Xb away from the ground plane. The different interferometer 
arms are separated from their neighbors by multiples of 9dZo, 
where 9d is the diffraction angle and zq is the longitudinal dis¬ 
tance to grating gl (in the case of the pillars and chopper cl) 
or grating g3 (in the case of chopper c2). a and R dimensions 
for the pillars and phase choppers are given in Table 1. 


grating in the ±x direction and observe the flux admitted 
through the third grating in order to determine the inter¬ 
ference pattern’s contrast and phase. We measure that 
transmitted atomic flux with a 100 ^m wide platinum 
wire Langmuir-Taylor detector [21]. 


In the rest of Section II we describe how we measure 
the atoms’ velocity distribution and polarizability. We 
measure vq, the atoms’ mean velocity, using phase chop¬ 
pers, which are charged wires parallel with the y axis held 
parallel to grounded planes, indicated in green in Fig. 1. 
We measure static polarizability a with a non-uniform 
electric field created by two oppositely charged cylindri¬ 
cal pillars parallel with the y axis and indicated in blue 
in Fig. 1. The pillars’ electric field shifts the interference 
fringe phase by an amount roughly proportional to ch/uq. 
Section IIA describes how the electric held geoemtry of 
both the phase choppers and the pillars causes a differ¬ 
ential phase shift. Section IIB describes how we use the 
phase choppers to measure the velocity distribution, and 
section IIC describes how we use the pillars to measure 
a. Section IID discusses how we apply our knowledge 
of the velocity distribution to analyze the polarizability 
data taken with the pillars. 


A. Phase shifts with cylindrical electrodes 


Both the pillars and the phase choppers are described 
by the geometry shown in Fig. 2, and create electric fields 
given by 


X — b X + b 

(x — 6)2 + (x + 6)2 + z^_ ^ 

_^_!_I5 (1) 

(a;-6)2 + z2 (x + b)2 + z2] 
where the effective line charge density 

exists a distance 6 = a^/T+~2Il/a away from the ground 
plane. The parameter a represents the distance between 
the ground plane and the closest cylinder edge, R repre¬ 
sents the pillars’ radius, and the directions x and z are 
shown in Fig. 2. 

When atoms enter an electric held, their potential en¬ 
ergy changes by f/stark = Since f/stark « -0.1 

/reV and F^kinetic ~ 1 eV for Cs in our experiment, we can 
use the WKB approximation along with the Residue The- 


E{x,z) = 


Sttco 
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orem to compute the total phase accumulated by an atom 
travelling through the field. We can also approximate 
that atoms travel parallel to the ground plane regardless 
of the angle at which they diffracted and their incident 
angle upon grating gl. Even though this approximation 
may be incorrect by up to 10“^ rad, such a discrepancy 
would only cause errors in the accumulated phases by 
factors of 10“®, which is insignificant for our experiment. 
Therefore, we represent the accumulated phase along one 
path for a component of an atomic de Broglie wave as 


<i)(w, x) 


1 

hv 



irelhv \b'^ — xl) 


( 3 ) 


where xt is the distance between the atom’s path and the 
ground plane. 

The atoms in our beam form many interferometers, 
but we only need to consider the four interferometers 
shown in Fig. 1. Other interferometers are insignificant to 
our analysis because they have some combination of low 
contrast and low flux. We label the four interferometers 
that we do consider with the index k = +2, fc = +1, 
k = —1, and k = —2. The differential phase shifts for 
the four interferometers are 


+ 26 'd 2 o) - ^{v,Xb + OdZo) 

= ^iv,Xb + OdZo) - ^{v,Xb) 

= ^{v,Xb) - ^iv,Xb - OdZo) 

^b) = Xb - OdZo) - 4>(u, Xb - 20dZo) (4) 

In the above equations, OdZo is the lateral separation be¬ 
tween classical paths in the interferometer, where Od = 
AdB /dg is the diffraction angle and zo the distance to the 
first grating (in the case of the pillars and chopper cl) or 
the third grating (in the case of chopper c2). 


B. Velocity measurement 

The atoms in the beam do not all have the same ve¬ 
locity, so the electric fields do not apply the same phase 
shifts to each diffracted atom. We observe the average 
phase and contrast of an ensemble of atoms with velocity 
distribution P{v). We model P{v) as a Gaussian distri¬ 
bution 


( 5 ) 

where vo is the mean velocity and the velocity ratio 
Vr = voluy is a measure of the distribution’s sharpness. 
It is worth noting that the velocity distribution for a su¬ 
personic atom beam is better described by a u^-weighted 
Gaussian distribution [13]. However, either distribution 
can be used in our analysis to parametrize the typical 
high-uo, high-Ur velocity distributions of our atom beam 
without changing our polarizability result by more than 




FIG. 3. (Color online) An example of a measurement of con¬ 
trast C vs. phase chopper frequency /c for a Cs beam. We 
fit a model to these data that has vo and Vr as fit parameters 
in order to measure the velocity distribution. 


0.008%. Since wq is the average velocity in a Gaussian 
but not in a u^-weighted Gaussian, we use Eqn. (5) to 
simplify our discussion of the error analysis. 

To measure vo and Vr, we use phase choppers [24, 25]. 
Each phase chopper is a charged wire about 1 mm away 
from a physical ground plane (see Table I for phase chop¬ 
per dimensions). Chopper cl is between the first two 
gratings and chopper c2 is a distance Zci,c 2 = 1269.3 ± 
0.25 mm downstream of chopper cl, between the last two 
gratings (see Fig. 1). The voltages on the choppers’ wires 
and the distances between the beam and the choppers’ 
ground planes are chosen such that chopper cl shifts the 
ensemble’s average phase by -|-7r and chopper c2 shifts it 
by -TT. 

When we pulse the choppers on and off at a frequency 
/c, an atom may receive a net phase shift of ±7r or 0 de¬ 
pending on its velocity and the time at which it passed 
through the first chopper. Holmgren et al. [24] gives an 
intuitive explanation of how we measure contrast C vs 
fc to determine vo and Vr- Fig. 3 shows an example of 
C vs fc data. Hromada et al. [25] later improved upon 
Holmgren et al.’s model of C vs fc by considering how 
the thickness and divergence of the beam causes some 
components of the atoms’ velocity distribution to not be 
detected. In the present work, we expanded our analysis 
to include the four interferometers shown in Fig. 1, per¬ 
formed a more in-depth error analysis, and added an ad¬ 
ditional calibration step to the measurement procedure, 
as we discuss next. 

Hromada et al. described how the thickness and diver¬ 
gence of the beam determines the likelihood for atoms of 
certain velocities to be detected [25]. The thickness and 
divergence is defined by the finite widths of the collimat¬ 
ing slits wi and W 2 ■ The finite width of the detector Wdet 
and the detector’s offset from the beamline in the x direc¬ 
tion Axdet also affect the probability of detecting atoms 
as a function of the atoms’ velocities and initial positions 
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in the apparatus. The phase and contrast we observe 
with our detector is that of an ensemble of atoms with 
different velocities, different incident positions on grating 
gl, and different incident angles on grating gl. 

Uncertainties in wi, W 2 , u>det, and Axdet are more sig¬ 
nificant for beams that are physically wider. In K beams, 
which have larger 64 of « 50 /irad and wider velocity 
distributions {vq ~ 2100 m/s, Vr « 14, and therefore 
ay « 150 m/s), more of the lower-velocity atoms in the 
distribution miss the detector. Therefore, uncertainties 
in the aforementioned quantities have a higher bearing 
on how we model the average velocity of detected atoms. 
Ignoring this component of the analysis would cause a 
systematic increase in measured vq by 0.5% and Vy by 
10% for a typical K beam. 

Modeling the four interferometers shown in Fig. 1, 
rather than only the k = ±1 interferometers, also im¬ 
proved our understanding of how likely it is for certain 
velocities to be detected. For K beams with wide veloc¬ 
ity distributions, we would report Vq too high by about 
0.5% and Vy too low by about 5% if we included only 
the k = ±1 interferometers. This is because, for such 
beams, a much higher proportion of atoms in the k = ±2 
interferometers miss the detector than in the fc = ±1 in¬ 
terferometers. Ignoring the k = ±2 interferometers has a 
significant effect on the model of the detected P{v) when 
the detected velocity distributions for the k = ±1 and 
k = ±2 interferometers are significantly different. Con¬ 
versely, for Cs and Rb beams, we found no significant 
difference in results between models because most of the 
atoms in all interferometers were detected regardless of 
velocity. 

Hromada et al. [25] also described how inequality 
between inter-grating distances Li and L 2 (see Fig. 1) 
causes systematic errors. When AL = L 2 — L 1 is nonzero, 
the interference fringes formed at the third grating be¬ 
come magnified or demagnified. We summarize this geo¬ 
metric magnification with the separation phase shift: 

Stt Z' k '\ 

A$sep,fc = ^^inc + 2 ^d,J AL (6) 

where k is the interferometer index (see Fig. 1) and O-mc 
is the incident angle on grating gl. To reduce systematic 
error in our vq and Vy measurements, we measure AL 
and set it equal to zero. 

Eqn. (6) implies that uncertainty in AL is more signif¬ 
icant for beams with larger 9d, such as K beams. Also, 
because A^sep,k has a component proportional to 0inc, 
uncertainty in AL is more significant for more divergent 
beams. As |AL| increases, uncertainties in wi, W 2 , w^et, 
and Accdet become more significant. Accordingly, we de¬ 
veloped a method to set AL = 0 to reduce those un¬ 
certainty contributions. Eqn. (6) implies that interfer¬ 
ometers on either side of the beamline receive opposite 
phase shifts. Therefore, by moving the detector in the 
±x direction, we observe linear changes in $ as a func¬ 
tion of Axdet with slope d$/da;det that is proportional 



-100 -50 0 50 100 


Ax^et UlTl) 



FIG. 4. (Color online) Two data plots showing how the phase 
>1> of the interference fringes changes linearly as a function of 
Aa;det and AL. In plot (a), the black dots represent coordi¬ 
nates at which data was acquired, and the colors represent 
contours inferred from the data. These figures show how $ 
vs Aa;det is a line with slope proportional to AL. 


to AL. Eig. 4 shows data that demonstrates this effect. 
We set AL to 0 ± 30 /rm by finding the AL for which 
d$/da;det = 0. 

Since we recalibrate AL every day, the 30 /rm uncer¬ 
tainty in AL represents a systematic error for one day’s 
measurements and a statistical error for many days’ mea¬ 
surements averaged together. That error will contribute 
toward the statistical uncertainty of a measurements. 
The same is true for Axdet = 0 ± 30 /rm, which also 
fluctuates from day to day as we set up the apparatus. 

If the interferometer grating bars are significantly non¬ 
vertical, it becomes necessary to consider the phase shift 
induced by the component of gravitational acceleration 
in the plane of the interferometer. That phase shift is 
given by 


Adlaccel — 


7rgsin(6lg)(Li -h L 2 Y 
2 dgV^ 


(7) 


where 9g is the tilt of the grating bars with respect to 
vertical [31, 32]. See [31] for an explanation of how we 
measured 9g. Our interferometer’s \9g\ never exceeded 
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□ Rb □ Cs 


^c1,c2 


= 30 ± 6 [im 
Wj = 40 ± 6 [im 
w^e, = 100 ± 3 |J.m 
= 1269.3 ± 0.25 mm 
3^1 = 986 ± 25 iJ.m 
ac 2 = 893 + 25 [im 
0 g = value ± 1.4 mrad 
total sys. error xIOOO 



0.0 0.1 0.2 


Svg/vo x 1000 5 v/V|. x 1000 


FIG. 5. (Color online) Systematic uncertainty budget for 
measurements of vo and Vr for our Cs, Rb, and K beams. 
The total systematic error in uo and Vr in turn contributes 
toward the total systematic uncertainty in a measurements 
(shown in Fig. 7). A nominal value for ffg is not listed be¬ 
cause ffg changed from —2.37±1.39 mrad to 1.73±0.59 mrad 
toward the end of the experiment. 


2.3 mrad. If we were to neglect this portion of the anal¬ 
ysis, we would report vq incorrectly by up to 0.015% and 
Vr incorrectly by up to 0.25%. 

The uncertainty budget for vq and Vr measurements 
is displayed in Fig. 5. The total statistical uncertainty 
in measured vg and Vr is roughly 10 times larger than 
the total systematic uncertainty after about 15 minutes 
of data acquisition with the phase choppers. Because vg 
and Vr drift over time, typically 3% over the course of 
several hours, we measure the velocity distribution twice 
every hour. 

Gaseous alkali atoms in an atomic beam nozzle have a 
probability of forming homonuclear dimers that depends 
on the gas pressure [33] and the diameter of the nozzle 
hole [34]. It is important for us to quantify the dimer 
mole fraction in our beam because the dimers’ spatially 
averaged (tensor) polarizabilities are approximately 1.75 
times the monomer polarizabilities [35]. In our nozzle, 
the vapor pressure of alkali atoms is on the order of 1 
torr at our typical running temperatures of 160° C for 
Cs, 220° C for Rb, and 350° C for K. According to data 
acquired by Gordon et al. (1971) [33] and Bergmann et 
al. (1978) [34], our alkali gas pressures of 1 torr should 
result in a dimer mole fraction well below 1%. Addition¬ 
ally, Holmgren et al. [11, 36] demonstrated how to place 
an upper limit on the dimer mole fraction by analyzing 
resolved diffraction patterns through a single nanograting 
and looking for peaks associated with dimer diffraction. 
In this work, we used very similar nozzle temperatures in 
our experiment as Holmgren et al. did in 2010. For all 
these reasons, we conclude that the dimer mole fraction 
in our beam must be less than 1%. Fig. 7 shows how a 
4% dimer mole fraction would lead to a significant (0.1%) 
error in measured polarizability. 


C. Polarizability measurement 


To measure the Cs, Rb, and K polarizability, we 
use two parallel, oppositely charged, 1/2-inch-diameter, 
stainless-steel pillars. The pillars are mounted to a sin¬ 
gle, rigid support structure so that a 3999.7 ± 1.0 ^m gap 
exists between them. A motor moves the support struc¬ 
ture in the ±a; direction, and a length gauge monitors 
the structure’s x position. The length gauge measures 
displacements of the structure with 30 nm accuracy. We 
begin a polarizability measurement with the assembly 
positioned such that the beam passes through the gap 
between the pillars near one of the edges. We take 25 sec 
of data with the electric field on and 25 sec with it off. We 
then move the pillars in nine 400 /rm increments so that 
the beam approaches the other edge of the gap, taking 50 
sec of data at each location. In doing so, we observe the 
phase shift A$ = ‘f’pinars.on — ‘href applied by the pillars 
as a function of Xb (see an example in Fig. 6). We then 
repeat this sequence, moving the pillars in the opposite 
direction in order to minimize possible systematic errors 
associated with travelling in a certain direction. When 
the electric field is off, we observe the reference phase 
‘href and reference contrast Gref given by 

1 1 "°° 

Grefe*'*’™‘ = Goe*‘^»-V / (8) 

2 ^ Jv=o 


The Sagnac phase, A<l>sag, is a phase shift caused by the 
Earth’s rotation and is described in [11, 37, 38]. Cg is 
the contrast that would be observed in the absence of 
A4>sag(u), and $0 is an arbitrary phase constant. When 
the field is on, we instead observe 


^‘^’pilla 


^pillars,on fc- 
1 


Goe 


i'I’n 


2 ^ 

k • 


P{v)e 


Av,Xb)+iAisa,g{v) 


dv (9) 


v—0 


We fit a model to A$ vs Xb, as shown in Fig. 6. The fit 
parameters of that model are the polarizability and the 
pillars position Xbo for which the phase shift is zero (i.e. 
the location of the virtual ground plane). 

In our earlier work, we used one pillar next to a 
grounded plate instead of two pillars forming a virtual 
ground plane [11]. We measured Xb by blocking the beam 
with the pillar. There were significant statistical errors 
of a few /im associated with this procedure, and a 1 /xm 
error would lead to a 0.1% error in polarizability. Our 
new pillars assembly greatly reduces those statistical er¬ 
rors. Measuring A4> vs Xb on both sides of the ground 
plane makes our typical 5 /xm uncertainties in Xbg add 
an insignificant amount of statistical uncertainty to the 
determined a. 

The systematic uncertainty budget for our polarizabil¬ 
ity measurements is shown in Fig. 7. In the next few 
paragraphs, we discuss how we measured some of the 
quantities in the error budget. 
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_c2 

0.1 

CO 




■o 

0.0 

'tfi 


0) 

-0.1 



Xb-Xbo(mm) 


FIG. 6. (Color online) An example of a measurement of phase 
shift vs X position of the pillars for a Rb beam. The two fit 
parameters used to fit the model to these data are polariz¬ 
ability aub and the pillars’ position at which the phase shift 
is null Xbo- 


TABLE II. A typical sequence of measurements during a day 
of data acquisition. The +x direction is arbitrarily chosen— 
the important aspect is that we spend an equal amount of 
time scanning the pillars in each direction so as to minimize 
possible systematic errors. This sequence of eight measure¬ 
ments is repeated once per hour for anywhere between 8 and 
36 hours. We end the data acquisition by repeating the first 
four measurements. 


Type of data acquired Duration 

contrast vs chopping freq. 7m 5s 

chopper cl phase 3m 45s 

chopper c2 phase 3m 45s 

contrast vs chopping freq. 7m 5s 


A4> vs pillars position {+x direction) 8m 45s 

A4> vs pillars position {—x direction) 8m 45s 

A4> vs pillars position {+x direction) 8m 45s 

A<I> vs pillars position {—x direction) 8m 45s 


■ K □ Rb □ Cs 


Wdet = 100 ± 3 pm 
dg = 99.90 ± 0.5 nm 
Vpiiiars: 8V/V = 0.0005 
apillars = 1999-9 ± 0.5 pm 
^pillars ~ 6350 + 0.5 pm^ 

Vq: Svg/Vo given in Fig. 5 
v,-: 5vr/V|. given in Fig. 5 
Zgi pillars = 833.50 + 0.25 mm 
dimer fraction < 0.01 
totai sys. error xIOOO 

0.0 0.4 0.8 1.2 
5a/a x1000 

FIG. 7. (Golor online) Systematic uncertainty budget for po¬ 
larizability measurements for our Cs, Rb, and K beams. The 
uncertainties in knowledge of vo and Vr are propagated for¬ 
ward from Fig. 5. Values of Vpmars ranged from 5 kV to 7 kV, 
always with 0.05% uncertainty. 


We reduced the uncertainty in tpuiars to 0.05% by in¬ 
dependently calibrating our voltage supplies. To measure 
V'piiiars, which ranged from 5 kV to 7 kV, we used a Fluke 
80K-40 high voltage probe. We measured the probe’s 
voltage divider constant, which itself depended on input 
voltage, using two Fluke 287 digital multimeters. 

We measured Zgi,piiiars to 1/4 mm accuracy. We placed 
rulers in the apparatus, after which three of us would read 
the rulers both live and from photographs. We repeated 
this process for many longitudinal positions of the rulers 
to further reduce statistical error in the measurement. Fi¬ 
nally, we compared the rulers we used with other rulers 
to verify that the ones we used were printed without sig¬ 
nificant systematic error. The value of Zgi.piiiars we use 


in this analysis is the average of all those measurements. 

We measured the width of the gap between the pil¬ 
lars, 2apiUars, to 1 p,m accuracy by repeatedly scanning 
the pillars assembly across the beam and recording the 
positions at which each pillar blocked half of the atom 
beam. To verify that OpiUars did not change over time, 
we repeated this procedure many times throughout the 
months during which we acquired our data. 

Our 9g was always close enough to zero such that we 
did not need to consider A<i>accei in our polarizability 
data analysis. We would only need to consider A$accei if 
\9g \ exceeded 23 mrad. We also find that uncertainties in 
wi and W 2 each do not correspond to more than 0.004% 
uncertainty in a. 


D. Determining the velocity distribution during 
polarizability measurements 

A typical sequence of measurements is shown in Table 
II. We measure the velocity distribution twice between 
every four scans of the pillars across the beam, and cal¬ 
ibrate the phase choppers between each pair of velocity 
measurements. 

We linearly interpolate between vq and Vj. measure¬ 
ments before and after each pillars scan to estimate those 
quantities at the time of that scan. Using cubic spline 
interpolation and Gaussian Process Regression to inter¬ 
polate between vq and Vj. measurements changes our re¬ 
ported polarizabilities by no more than 0.001%, which is 
small compared to our other uncertainties. 


III. RESULTS AND DISCUSSION 

Table VII shows our measurement results for the K, 
Rb, and Cs atomic polarizabilities, Ok, <aRb, and acs- 


























TABLE III. Direct measurements of Cs, Rb, and K static, 
ground-state polarizabilities. 


Atom 

a(stat.)(sys.) (A^) 

Gs 

59.39(3)(9) 

Rb 

47.39(3)(8) 

K 

42.93(2)(7) 


TABLE IV. Measured ratios of Cs, Rb, and K static, ground- 
state polarizabilities. The systematic uncertainties in each 
ratio, which arise from the fact that the systematic errors 
in different measurements are not perfectly correlated, are 
negligible compared to the statistical uncertainties. 


Ratio 

Value (st at.) 

Sys. Err. 

acs/ttK 

1.3834(9) 

3•10-” 

aCs/aRb 

1.2532(10) 

7-10-^ 

OCKh/oCK 

1.1040(9) 

2•10"^ 


The tabulated statistical uncertainties are the standard 
error of the mean for each result. To get this statistical 
precision, we acquired over 90 hours of data, including 
150 data sets similar to Fig. 6 and 60 data sets similar to 
Fig. 3. The total systematic uncertainty for each mea- 
snrement is also stated in Table VII, and a breakdown 
of the systematic uncertainty budget is summarized in 
Fig. 7. While the statistical uncertainties are typically 
0.05% for our measurements of polarizabilities, the sys¬ 
tematic uncertainties are three to four times larger, and 
cause a total uncertainty of typically 0.16% for each mea¬ 
surement. 

We report the ratios of polarizabilities q;cs/q;kj 
aCs/cTRb, and anh/oiK in Table IV. These ratios have 
uncertainties smaller than 0.08% because we used the 
same apparatus for each direct measurement. For many 
of the sources of systematic uncertainty summarized in 
Fig. 7, an error in one of those quantities would scale each 
direct polarizability measurement by the same amount. 
These correlated uncertainties, such as electrode geom¬ 
etry or grating pitch, do not contribute significantly to 
systematic errors in our measured polarizability ratios. 
However, uncertainties in ici, W 2 , Wdet, and AL affect our 
Q!k, OfRb, and acs measurements differently and therefore 
contribute a small amount to the final systematic uncer¬ 
tainties in the ratios. Even so, the ratios’ systematic 
errors are much smaller than the statistical uncertain¬ 
ties. We discuss the value of high-precision ratios for 
testing atomic theories in Section III A, and we discuss 
the possibility of using such ratios to improve individual 
measurements of polarizability in Section IV). 


A. Comparisons with other experimental and 
theoretical polarizabilities 

Fig. 8 and Table VIII compare our polarizability mea¬ 
surements with ab initio calculations, semi-empirical cal¬ 
culations, and experimental measurements subsequent to 
and including Molof et o7’s and Hall et oL’s 1974 mea¬ 
surements [4, 7]. First, we will discuss the comparison 
to previous direct measurements. Our K and Rb polar¬ 
izability measurements have 3 times smaller uncertainty 
than our group’s previously published direct measure¬ 
ments of ax and aub [H]) and 10 times smaller uncer¬ 
tainty than the only other direct measurements of ax 
and axb, which were made using the E-H gradient bal¬ 
ance technique [7] and the E gradient deflection technique 
[4]. We emphasize that our new measurements are inde¬ 
pendent of the results in [11] because although we used 
the same atom interferometer machine, we used a dif¬ 
ferent material nanograting gl, different electrodes with 
different geometry, a different atom beam velocity mea¬ 
surement technique, a different atom beam source noz¬ 
zle, and a detector with a different width. Hence, the 
fact that our new and more precise measurements are 
consistent with the measurements in [4, 7, 11] should be 
regarded as an independent validation of each of these 
previous results. 

There is one other direct measurement of acs with un¬ 
certainty similar to (and slightly smaller than) ours. Our 
direct acs measurement is 11 times more precise than 
the result reported using the E-H balance technique [7], 
but 1.1 times less precise than Amini and Gould’s 2003 
measurement [8] that was made using an atomic foun¬ 
tain apparatus. To our knowledge, Amini and Gould’s 
work is the only polarizability measurement to date that 
has been accomplished using an atomic fountain, and 
it produced a remarkable improvement in precision by 
a factor of 15 as compared to the only previous direct 
measurements of acs [4, 7]. Furthermore, acs measure¬ 
ments can test some of the atomic structure theory that 
is used to interpret atomic parity non-conservation ex¬ 
periments as a way of constraining physics beyond the 
standard model [15-18]. Thus, it is particularly impor¬ 
tant to validate this acs result in [8]. We find that our 
acs measurement is consistent with Amini and Gould’s. 
Our result acs = 59.39(9)A deviates from their result 
of acs = 59.42(8)A^ by 0.03 A^, which is insignificant. 
Gomparing our atom interferometer result with their 
fountain result serves as a cross-check for both methods. 
Both measurements also agree with acs values inferred 
from the atomic structure calculations by Derevianko and 
Porsev (2002) [17] and Derevianko et al. (1999) [50] for 
PNG analysis. 

Most theoretical predictions for ax: ^xb: and acs de¬ 
viate from each other and from our measurements signifi¬ 
cantly. Out of 28 sets of theoretical predictions shown in 
Fig. 8, only ten sets of predictions [17, 50, 55-62] are 
consistent with our results within 3a (where cr is the 
standard deviation of our measurement). Furthermore, 
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FIG. 8. (Color online) Our direct measurements (a) and measured ratios (b) compared with other measurements, ab initio 
calculations, and semi-empirical calculations [4, 7, 8, 11, 17, 39-63]. The references are represented on the x-axis by the 
first three letters of the first author’s last name followed by the year of publication. For the semi-empirical calculations: 
Reference Fue82 used semi-empirical pseudopotentials [42], Pat97 used experimentally-determined energy levels [48], Der99 
used experimentally-determined electric dipole transition matrix elements [50], and Mit03 and Yonl4 used semi-empirical core 
polarization potentials [52, 63]. Values in this plot from 1999 and later are also reported in Table VIII at the end of this paper. 


the semi-empirical ax, ofRb, and acs values calculated 
in 1999 by Derevianko et al. [50] are the only predic¬ 
tions that match all three of our own ctk, ctRb, and acs 
measurements to within Scr. These predictions [50] were 
made using measured lifetimes and energies, and Dere¬ 
vianko and Porsev’s later acs prediction [17] was made 
using a measured van der Waals Cq coefficient. This 
is an important point because there are now additional 
data on lifetimes, van der Waals Cq measurements, and 
line strength ratios that can inform new semi-empirical 
predictions for polarizabilities that we discuss in Section 
IIIB and Fig. 9. 

Fig. 8 (b) compares our measurements of atomic polar¬ 
izability ratios to other theoretical, semi-empirical, and 
experimental reports for these ratios. The values we 
measured for acs/aK, ctCs/aRb, and Q;R,b/<aK are consis¬ 
tent with all of the previous experimental measurements 
of these ratios, given the larger uncertainties associated 
with previous measurements. Comparing theoretical pre¬ 
dictions to our measured polarizability ratios serves as a 
different way to test the theoretical predictions. Since 
the fractional uncertainties on our measured ratios are 
smaller than those of our absolute measurements, our 


ratios serve as a more precise test for theoretical works 
that predict a values for multiple alkali atoms. 


B. Comparisons with polarizabilities derived from 
other quantities 

Static polarizabilities can be related to electric dipole 
transition matrix elements, state lifetimes, oscillator 
strengths, and van der Waals coefficients. We will de¬ 
scribe those relations and compare our a measurements 
to (X values derived from recent calculations and high- 
precision measurements of those quantities. Those com¬ 
parisons are shown in Fig. 9 and Table VIII. 

The polarizability (in volume units) of an atom in state 
i can be written in terms of Einstein A coefficients 
as 


a 


_ ^ Ajj gj 

* - 9 ,„4. 




( 10 ) 


where is the transition frequency between states i 
and k, and = 2J„ -|- 1 is the degeneracy factor for 
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FIG. 9. (Color online) Comparisons of our lab’s polarizabil¬ 
ity measurements (this work as well as [11]) and Amini and 
Gould’s acs measurement [8] to polarizabilities derived from 
measured lifetimes and lifetime ratios, [64-73], lifetimes in¬ 
ferred from photo-association data [74, 75], theoretical 
values [76], and van der Waals Ce measurements ]17, 77, 78]. 
Values in this plot from 1999 and later are also reported in 
Table VIII at the end of this paper. 


TABLE V. We use the following residual polarizabilities Or 
and matrix element ratios R = 713 / 2 /^ 1/2 • The sources for 
each quantity are cited next to the values in the table. 


Atom 

ar (A^) 

R 


Gs 

2.481(16) [17] 

1.9809(9) [73] 

Rb 

1.562(89) [79] 

1.996(4) 

[65] 

K 

0.925(45) [79] 

2.000(4) 

[80] 


state n. In our case, state i is the ground state. The 
residual polarizability ttr includes terms not explicitly in¬ 
cluded in the sum, the polarizability of the core electrons, 
and a correction accounting for correlations between core 
and valence electrons as described in several references 
[17, 50, 79]. We will explicitly sum over the principal 
transitions from nsi /2 to npi /2 and np 3 / 2 , where n = 6 
for Cs, n = 5 for Rb, and n = 4 for K, and we will 
include the other transitions in ar- We abbreviate the 
lifetimes associated with the principal transitions to ri /2 
and T 3 / 2 . In our calculations, we use the transition wave¬ 
lengths uJik from references [81-84] and the ar values indi¬ 
cated in Table V. Fig. 9 shows polarizabilities calculated 
using measurements of T 1/2 and T 3/2 [64-66, 68-70, 72]. 
Fig. 9 also shows ocs calculated from values of Ti/ 2 _Rb 
and T 3 / 2 _Rb inferred in 2002 by Gutterres et al. from 
photo-association data taken in 2000 by Gabbanini et al. 
[74, 75]. 


We can use Patterson et al.’s 2015 measurement of 
'’ 3 / 2 ,Cs [85] along with a measurement of the ratio of prin¬ 
cipal transition matrix elements to report acs- Rafac and 
Tanner measured the ratio of Gs electric dipole transition 
matrix elements [73] 


Rcs 


( 6si/2 || i/|l 6P3 / 2 ) 
( 6si/2 || T >|1 6P1 / 2 ) 


which is related to the ratio of lifetimes 

Ti /2 _ Rcs ( 1^3/2 
T3/2 2 \Wi/2 


( 11 ) 


( 12 ) 


We can also report a polarizability using i?cs [73] 
in conjunction with Porsev et al.’s 2010 calculation of 
D ^/2 Cs ~ 20.334 (in atomic units) [76]. We can write ai 
in terms of the electric dipole transition matrix elements 
as 


CXi — 


12776000 


E 


1^11 J) 


Huji 


+ ar 


(13) 


where oq is the Bohr radius. As before, we only explic¬ 
itly consider the nsi /2 — npi /2 and nsi /2 — 0 ^ 3/2 matrix 
elements, where nsi /2 is the ground state. We abbre¬ 
viate the matrix elements associated with the principal 
transitions to D 1/2 and I? 3 / 2 . 


In 2002, Derevianko and Porsev demonstrated a 
method for obtaining values of and cs fro™ 

Gs van der Waals Cg coefficients [17] and Rcs [73]. Fig. 9 
includes a values derived using experimental Gs Cq mea¬ 
surements in conjunction with that method [77, 78]. 
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TABLE VI. Matrix elements, lifetimes, oscillator strengths, 
line strengths, and van der Waals Ce coefficients calculated 
from our polarizability measurements. We used R values from 
Table V. The matrix elements, line strengths, and Cq coef¬ 
ficient are expressed in atomic units, while the lifetimes are 
expressed in SI units. 5a, Sr, and Sa,. represent the uncer¬ 
tainties in the values due to uncertainty in a, R, and a^, 
respectively, dtot is the total uncertainty in the value. (-) 
represents a nonzero uncertainty of less than 0.5. 


Quantity Atom Value 

6a 

6r 

6a,. 

6tot 

Di/2 

Cs 

4.508 

(4) 

(1) 

(1) 

(4) 

Rb 

4.239 

(4) 

(3) 

(4) 

(6) 


K 

4.101 

(3) 

(3) 

(2) 

(5) 

D 3/2 

Cs 

6.345 

(5) 

(-) 

(1) 

(5) 

Rb 

5.989 

(5) 

(2) 

(6) 

(8) 


K 

5.800 

(5) 

(2) 

(3) 

(6) 

Ti /2 (ns 

0 Cs 

34.77 

(5) 

(1) 

(1) 

(6) 


Rb 

27.60 

(5) 

(4) 

(5) 

(8) 


K 

26.81 

(4) 

(4) 

(3) 

(6) 

T 3/2 (ns 

0 Cs 

30.37 

(5) 

(-) 

(1) 

(5) 


Rb 

26.14 

(5) 

(2) 

(5) 

(7) 


K 

26.45 

(4) 

(2) 

(3) 

(6) 

fl/2 

Cs 

0.3450 

(5) 

(1) 

(1) 

(6) 

Rb 

0.3433 

(6) 

(5) 

(7) 

(10) 


K 

0.3317 

(6) 

(4) 

(4) 

(8) 

/3/2 

Cs 

0.7174 

(11) 

(1) 

(2) 

(12) 


Rb 

0.6982 

(12) 

(5) 

(14) 

(19) 


K 

0.6665 

(11) 

(4) 

(7) 

(14) 

'S'l/2 

Cs 

20.32 

(3) 

(1) 

(1) 

(3) 


Rb 

17.97 

(3) 

(3) 

(3) 

(5) 


K 

16.82 

(3) 

(2) 

(2) 

(4) 

S 3/2 

Cs 

40.26 

(6) 

(1) 

(1) 

(6) 

Rb 

35.87 

(6) 

(3) 

(7) 

(10) 


K 

33.64 

(6) 

(2) 

(4) 

(7) 

Cq 

Cs 

6879 

(20) 

(-) 

(7) 

(21) 


Rb 

4719 

(15) 

(-) 

(26) 

(30) 


K 

3884 

(13) 

(-) 

(14) 

(19) 


C. Other atomic properties derived from our 
polarizability measurements 


Finally, we use our polarizability measurements to re¬ 
port matrix elements, lifetimes, oscillator strengths, line 
strengths, and van der Waals Cq coefficients. In these 
calculations, we use residual polarizabilities and ma¬ 
trix element ratios R from Table V. To report matrix 
elements and lifetimes, we use Eqn. (13) and Eqn. (11). 
at is given in terms of oscillator strengths fij as 


a, = 


E- 

3¥=^ 


(14) 


TABLE VII. Excited state polarizabilities where n = 

6 for Cs, n = 5 for Rb, and n = 4 for K. The values were 
calculated using our measurements and mea¬ 

surements [86, 87]. 


Atom 


Cs 

196.81(9) 

Rb 

120.33(8) 

K 

89.92(7) 


where m is the electron mass, is also given in terms 
of line strengths Sji as 


E 


S'. 


3^ 


dTTeoh ^ giOJij 


(15) 


Cq can be expressed in terms of dynamic polarizability 
as 

o roc, 

Cq = — a{iuj)'^duj (16) 

Jo 


Derevianko et aids 2010 work tabulates values of a{iuj) 
for Cs, Rb, and K atoms among others [88]. To derive 
Cq values from our q;( 0) measurements, we modify Dere¬ 
vianko et aids values of a{iuj) to get 


0^new(^^) — d(fcc) [cTthis work(O) 


am 


ap{ioj) 

ap{ 0 ) 


(17) 


where a{iuj) and d(0) refer to values tabulated by Dere¬ 
vianko et al. In the above equation (17), ap{iuj) is the 
contribution to a(iuj) by the principle transitions. The 
ratio ap(iuj)/ap{ 0 ) is given by 


ap(iuj) 

ap{0) 


-^1/2 

•*■’3/2 


^1/2 


I 

1/2 


3/2 




R 


^ 1/2 

"^ 3/2 


(18) 


Predictions of the parity-non-conserving amplitude, 
EpNC) in Cs depends heavily on Di/ 2 . We note that 
our Cs Di /2 value is consistent with the theoretical Cs 
Di /2 calculated by Porsev et al. in 2010 for the purpose 
of interpreting PNC data as a test of the standard model 
[76]. 

Finally, we use our measurements together with recent 
measurements of Cs, Rb, and K — o:nsi /2 [86, 87] 

to report excited state polarizabilities cxnp^^^ with better 
than 0.08% uncertainty. These results are shown in Table 
?? and serve as benchmark tests for calculations of dipole 
transition matrix elements for p — d transitions. 
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TABLE VIII. Comparisons of directly measured, semi-empirical, and ab initio polarizabilities from 1999 and later. These values 
are also plotted in Fig. 8 and Fig. 9. Lines with Or in the Method column were calculated using residual polarizabilities Or 
in Table V [17, 79]. Ratios of polarizabilities Q;Rb/«K = 1.097(5) from [11] and acs/oRb = 1.2532(10), acs/oK = 1.3834(9), 
ttRb/oK = 1.1040(9)from this work have been reported with smaller fractional uncertainties than the direct measurements 
listed in this table. For the ab initio method acronyms: CCSD indicates the coupled-cluster approach with single and double 
excitations of the wavefunctions. The prefix R stands for “relativistic” and the prefix L stands for “linearized”. The suffix T 
indicates that the authors used perturbative triple excitations of the wavefunctions. See ]1] for explanations of these various 
theoretical methods. 


Reference(s) 

Method 

acs (A ) 

aRb (A ) 

ttK (A"') 

Raf99 [17, 69] 

U/ 2 ) 7 - 3/2 nieas. -b 

58.97(22) 



Der99 [50] 

ab initio, RLGGSD 

59.50 

46.89 

42.84 

Der99 [50] 

semi-empirical 

59.26(28) 

47.21(9) 

43.00(12) 

LeoOO [17, 73, 77] 

Cq meas. + Rcs meas. + thry -|- 

59.49(25) 



Gut02 [74, 75, 79] 

Ti/ 2 , T 3/2 from PA data -I- ar 


47.18(22) 


Der02 [17] 

semi-empirical 

59.35(12) 



Mag02 [51] 

ab initio 

64.31 

49.64 

44.75 

Ami03 [8] 

direct acs meas. 

59.42(8) 



Mit03 [52] 

semi-empirical 


46.78 

42.97 

Ghi04 [17, 73, 78] 

Cq meas. -|- Rcs meas. -1- thry + ar 

59.35(22) 



Saf04 [53] 

ab initio 

59.00(13) 



Lim05 [54] 

ab initio, RGGSDT 

58.55 

47.29 

43.08 

Fal06 [70, 79] 

''‘ 1 / 2 , 7 - 3/2 meas. -b ar 



43.02(8) 

Bou07 [17, 72] 

Til’ll T-iji meas. -b ar 

59.31(82) 



Isk07 [55] 

ab initio, RLGGSDT 

59.04(10) 



Saf08 [56] 

ab initio, RLGGSDT 



42.87 

HollO [11] 

direct aRb and ctk meas. 


47.24(44) 

43.06(36) 

HollO [9, 11] 

ratio calibrated with aNa 


47.24(21) 

43.06(21) 

PorlO [17, 73, 76] 

ab initio -b i?cs meas. -b ar 

59.42(3) 



Safll [57] 

ab initio, RGGSD 


47.72(59) 


Nanl2 [58] 

ab initio, RGGSDT 



43.05(15) 

Borl3 [59] 

ab initio, RGGSDT 

59.13 



Safl3 [60] 

ab initio, RLGGSDT 



43.03(9) 

Sahl3 [61] 

ab initio, RGGSDT 



42.94(9) 

Jail3 [62] 

semi-empirical 



42.98 

Yonl4 [63] 

semi-empirical 

58.72 

47.07 

42.94 

Patl5 [17, 73, 85] 

T 3/2 meas. -b i?cs meas. -b ar 

59.21(13) 



This work 

direct meas. 

59.39(9) 

47.39(8) 

42.93(7) 


IV. OUTLOOK 

We are currently exploring ways to measure the po¬ 
larizability of Li and metastable He, the polarizabilities 
of which can be accurately calculated. By measuring 
ctCs/ctHe* or acs/oLi, we could report acs with preci¬ 
sion comparable to that of the ratios reported here for 
the benefit of PNC research. Such a measurement would 
also act as a calibration of the measurements presented in 
this work, because it would be independent of systematic 
errors that may affect our direct measurements. 


We are also exploring electron-impact ionization 
schemes for atom detection, which would allow us to de¬ 
tect most atoms and molecules. Our Langmuir-Taylor 
detector only allows us to detect alkali metals and some 
alkaline-Earth metals [21]. Installing an electron-impact 
ionization detector would allow us to broaden the scope 
of atom interferometry as a precision measurement tool. 
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